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ABSTRACT 

We report HST/NICMOS coronagraphic images of the HD 15115 circumstellar disk at 1.1 /im. We find a 
similar morphology to that seen in the visible and at H band-an edge-on disk that is asymmetric in surface 
brightness. Several aspects of the 1.1 /im data are different, highlighting the need for multi-wavelength images 
of each circumstellar disk. We find a flattening to the western surface brightness profile at 1.1 /im interior to 2" 
(90 AU) and a warp in the western half of the disk. We measure the surface brightness profiles of the two disk 
lobes and create a measure of the dust scattering efficiency between 0.55-1. 65/im at 1", 2", and 3". At 2" the 
western lobe has a neutral spectrum up to 1.1 /im and a strong absorption or blue spectrum > 1 . 1 /tm, while a 
blue trend is seen in the eastern lobe. At 1" the disk has a red Fl lOW-H color in both lobes. 
Subject headings: stars :individual(HD 151 15)-circumstellar matter 
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1. INTRODUCTION 

Raw material from the interstellar medium is processed 
through circumstellar disks into many different types of plan- 
etary bodies. Debris disks are a useful environment to study 
the composition of planetesimals that form around stars dif- 
ferent from our own Sun. Multi-wavelength scattered light 
observations of disks provide one avenue for determining the 
compositon of dust caused by collisions of planetesimals. 

HD 15115 possesses one of a growing number of spatially 
resolved debris disks amenable to multiwavelength observa- 
tions. HD 15115 is an F2 star at a distance of 45 ±1 pc 
jvan LeeuwenI |2007|) a nd was observed to have an IR ex- 
cess (I Silverstonell2000h . Unresolved emission from its disk 
was detected at 60, 100, and 850 iim, implying a ring at 
^35 AU with a temperature of 62 K (Zuckerman & Sona 
I2004t iDecinetal. 2003; Williams & Andrews 2006). Re- 
cently, iKala^^l] (120071 hereafter KFG07) discovered cir- 
cumstellar emission from a combination of Hubble Space 
Telescope/Advanced Camera for Surveys (HST/ACS) and 
Keck near-IR AO coronagraphic imaging. Surface bright- 
ness (SB) profiles were reported for V and H bands. Fur- 
thermore, they observed an extreme size and brightness dif- 
ference between the two lobes of the disk and shape asymme- 
tries around the disk midplane similar to that seen around f3 
Pictoris and indicative of a second disk at a different inclina- 
tion (Kalas & Jewitt 1995; Golimowski et al. 2006). We have 
imaged HD 15115 at 1.1 /im with HST; we present our obser- 
vations in ^ and analyze the results in ^ We discuss the 
implications of our work in ^ 

2. OBSERVATIONS 

We observed HD 15115 with the Hubble Space Tele- 
scope's Near-Infrared and Multi-Object Spectrometer (NIC- 
MOS) ([Thompson et al.ir998i) using the coronagraph in Cam- 
era 2. We used the Fl lOW filter with a central wavelength of 
l.l/im(0.8/im-1.4/im). Observations were taken on 2006 Oc- 
tober 20 for HD 15115 and HD 16647, the PSF reference. 

Accurate disk photometry requires a correct scaling ratio 
between the target star and its PSF reference. Both stars 
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would saturate the detector in the shortest exposure. There- 
fore, we estimated the Fl lOW flux density for both stars with 
synthetic photometry using the CALCPHOT task of the S YN- 
PHOT' package in conjunction with Kurucz models tabu- 
lated in the SYNPHOT library. We found a best fit match 
to HD 151 15's luminosity and effective temperature based on 
its parallax from the Hipparcos satellite, its Tycho-2 B and V 
photometry, and 2MASS /, H, and photometry. The model 
assumed a solar metallicity and the best fit was for a logg=5.0, 
Te//=6750 K star, with a luminosity of 3.11L0. We fit so- 
lar metallicity models to HD 16647, finding a best fit for a 
logg=5.0 Tf//=6750 K star and a luminosity of 4.0 IL©. The 
resulting SYNPHOT scali ng in F llOW is 0.600. Based on 
previous experience, (De bes et al.l l2008, for example), SYN- 
PHOT is accurate to ^3%, which we will adopt as our a priori 
estimate of the scaling uncertainty. Coronagraphic images of 
the disk were taken at two spacecraft orientations. Th images 
were subtr acted, registered, an d combined following the pro- 
cedures of iDebes et al.l (l2008l) . The individual subtracted and 
combined images are shown in Figure[T] 

3. ANALYSIS 

Figure [T] shows that the disk is detected at the same posi- 
tion angle (PA) as KFG07, despite the presence of one strong 
residual streak close to the position of the eastern lobe of the 
disk (Feature A in Figure[T]i. Perpendicular to the disk, strong 
residuals also exist. We measured the PA, full-width at half- 
maximum (FWHM) of the disk midplane, and SB profiles as 
a function of radius from the central star. 

3.1. Disk Geometry 

We measured the PA and FWHM of the western lobe of the 
disk from 0."68 (30 AU) to 3."69 (165 AU). For the eastern 
lobe we took measurements from 0"68 to 3"0 (135 AU). In 
our analysis we rotated the final disk images so that the disk 
midplane was oriented along an image row. We took 3x13 
(0"23 xO."98) pixel cuts of the disk centered on the brightest 
pixel of the disk and summed in a direction perpendicular to 
the midplane direction. We then fit a Gaussian curve to the 
disk emission to measure the location of the midplane (and 
thus the local PA) as weU as the FWHM of the disk. We 
took an average of the values for both orientations and es- 
timated the uncertainty as being half the difference between 

' http://www.stsci.edu/resources/software_hardware/stsdas 
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Fig. 1 . — (top) False color images of the observations of HD 15115 in the 
FllOW filter at two spacecraft orientations. The images are logarithmically 
scaled. (bottom)False color, final combined image of the HD 15115 disk, 
showing the contaminating feature A. The dashed line represents the nominal 
position angle of the disk at 278.9°, while the asterisk represents the position 
of the star. We have masked out separations in the disk within Of'68, where 
PSF subtraction residuals dominate. Noticeable are a warp on the western 
lobe, and the asymmetries between the western and eastern lobes. 



the two measurements. Figure |2] shows the PA and FWHM 
of the western lobe of the disk as a function of radius. Be- 
yond l."7, we find a PA for the disk of 278.9° ±0.2°, which is 
consistent with KFGOV's value of 278.5° ±0.5°. However, in- 
terior to 1 ."7, we note that the disk PA increases with decreas- 
ing distance from the star to a maximum PA of 290° ± 1 ° at 
0."68. The changing PA of the midplane is seen at both orien- 
tations. Furthermore, we find at 2", a FWHM of 0"29±0.04 
which is consistent with the FWHM of 0" 19±0. 1 reported by 
KFG07. The median FWHM of the disk from 0"68-3."69 is 
0."26±0"07. The FWHM at 1 . 1 fim is broadened by the NIC- 
MOS Fl low PSF, which has a FWHM of 0"1, implying an 
intrinsic FWHM for the disk midplane of '-0"24 or 1 1 AU. 

Measurement of the eastern side of the disk is complicated 
by the contamination of Feature A from PSF subtraction. To 
measure the FWHM and PA, we masked out pixels heavily 
contaminated by this feature. We again took the average of 
the two sets of measurements and estimated the uncertainties 
as being half the difference between the two measurements. 
We find no evidence of a warp in the eastern lobe, which has 
a median PA of 96.6° ± 1.3°. To compare to the PA of the 
Western side we add 180° to get 276.6°, consistent with no 
significant bowing i n the large scale structure of the disk as is 
seen for HD 61005 (Hines et al.ll2007l) . The median value of 
the eastern FWHM is 0."29±0"06. 

3.2. Surface Brightness Profiles 

Figure [3] shows a comparison of the 1 . 1 fim SB profiles in 
Am/71 loiv arcsec"^ relative to the central star (m/riiovv=6.12) 
for the eastern and western lobes of the disk by taking 
0"23x0."23 apertures centered on the midplane of the disk. 
Between 0.7-1.8", the disk SB of the western lobe falls as 
r-i-4±o i, while at radii >1.8", the SB drops as r-3-56±0 06 j^ie 
eastern lobe drops steadily as ^-2 i4±o.o6 fj-Q^ 0."7-3". 

3.3. Scattering Efficiency vs. Wavelength 

Combining our data with the reported V and H SB data, 
a measure of the scattering efficiency of the dust around 
HD 15115 can be constructed. The measured SB at any ra- 
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Fig. 2. — Measurements of the HD 15115 disk midplane PA (top) and 
FWHM (bottom) as a function of distance for the two lobes of the disk. 
Squares represent the western lobe, while triangles represent the eastern lobe. 
The sharp rise in PA for the western lobe interior If 7 is due to the warp, and 
the Eastern lobe's PA beyond 1 " 8 is marginally different than the Western 
lobe. 



dius in the disk can be self-consistently modeled for com- 
position and grain size distribution provided a knowledge of 
the density distribution of the dust is known (see details in 
iDebes et aLl i2008) . The SB of a dust disk at a particular ra- 
dius is oc F^,(f>(9)Qscaj^, where 0(6*) is the phase function of 
the dust at a particular scattering angle, Qsca is the scattering 
efficiency of the dust, and is the ratio of stellar luminosity 
to the IR luminosity of the dust distribution. This is difficult 
for HD 15115 because the density distribution of the dust is 
not known through, for example, spatially resolved thermal 
emission images. However, multi-wavelength imaging can 
provide useful compositional information even in the absence 
of spatially resolved thermal emission constraints. In the case 
of HD 15115, three images at different wavelengths will be 
insufficient to provide a definitive picture, but can constrain 
certain compositions. 

In order to compare results across different wavelengths, 
instruments, and telescopes, it is imperative that a careful 
and consistent treatment of SB measurements is implemented. 
Observations of disks at different wavelengths can be affected 
by the respective PSFs of telescopes with different diffrac- 
tion limits and Strehl ratios, but no standard technique is 
applied to the study of circumstellar disks. Intuitively, one 
would expect minor corrections as a function of wavelength 
provided that the disk structure did not significantly change 
between wavelengths and that the disk itself was mostly re- 
solved. KFG07 estimated aperture corrections using point 
source PSF s. Other methods have included deconvolution 
of the data jGolimowski et alj|2006l) . or using apertures that 
capture the majority of the observed emission from the disk 
at each waveleng th (Schn eider et al.l [19991 : [Weinberger et alJ 
I2OO2I: iDebes et al. 2008). We modeled the effect of the dif- 
ferent PSFs of NICMOS, ACS, and Keck H-band adaptive 
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Fig. 3. — The surface brightness profiles of the eastern and western lobes 
of the disk in FllOW magnitudes per square arcsecond for a 0f'23x0f'23 
aperture. No connections were applied. The symbols are the same as in Figure 



optics on the measured surface brightness for the 0"23 x0"23 
aperture. 

Our approach to calculating the corrections derives from 
the reliable theoretical PSFs generated by TinyTim^ for HST 
as well as using a Keck H PSF corresponding to the October 
observations of HD 15115 (M. Fitzgerald, personal communi- 
cation). We took the PSFs for each instrument and convolved 
them with model disks constructed from the measured power- 
law slopes and disk FWHMs for each observation. We then 
measured the SB as a function of distance for the original and 
convolved models with the 0."23x0"23 aperture in FllOW 
and a 0"25 x0"25 aperture for V and H and took the ratio as 
a single valued multiplicative correction to the observed SB 
profiles. We verified that the slope of the surface brightness 
profiles did not change due to convolution. To estimate the 
uncertainty in this correction, we used the standard deviation 
of the measurements between 2"-3", the region where all the 
observations overlap. For the ACS coronagraphic PSF, we 
found a correction of 1.23 ±0.02, for the H band we found a 
correction of 1.56±0.06, and for FllOW we found a correc- 
tion of 1.36±0.04. 

To compare the NICMOS and other measurements, we 
must make a correction for the slightly different photomet- 
ric aperture sizes of 0."23x0."23 and 0"25x0"25 respec- 
tively. To do this we oversampled our model NICMOS im- 
ages by a factor of 3 and calculated the predicted difference 
in SB measurements between our 0."23x0."23 aperture and a 
0"25 x0."25 one. This multiplicative correction changes with 
distance from the star, and ranges from 1.2 to 0.8. 

In Figure |4] we plot SB/Fj^ ^, as a function of wavelength 
for HD 15115's two disk lobes at 1", 2", and 3". This ob- 
served quantity is oc Qsca- Our 1.1 /im data provides a third 
point in the measure of HD 15 1 15's scattering efficency when 

^ http ://w w w. stsci . edu/software/tinytim/ 
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Fig. 4. — Comparison of the scattering efficiency (or color) vs. wavelength 
of HD 15115's disk at 1"(45 AU), 2"(90 AU), and 3"(135 AU). The sym- 
bols are the same as in Figure |2] The disk dramatically changes color over 
wavelength and distance from the central star. 



combined with the 0.55/xm (V) and 1.65/xm (H) measure- 
ments in KFG07. We derive the uncertainties in the scatter- 
ing efficiency for 0.55/xm and 1.65/im based on Figure 3 in 
KPG07. Where no error bars were given for the 1.65/im data 
in KFG07, we assumed an uncertainty of 20%. In the western 
lobe at 2", we find that the disk's scattering efficiency is neu- 
tral out to 1 . 1 /im but then drops such that the spectrum is blue 
at 1.65/im. The spectrum is the same from 0.55-1.65/tm at 
3". At 2", the eastern lobe is strongly blue and at 3" the blue 
trend continues out to 1.1 /im. KFG07 did not report SB for 
the H data at 3", but if the spectrum is the same as the inner 
part of the disk, one would expect the blue trend to continue 
to longer wavelengths. 

Finally, interior to about 1 "8, we can get a rough sense of 
the near-lR colors of the disk in the east and west. We took 
KFG07's measurements of the two lobes in H at l"and com- 
pared them to our data. In both cases the disk is red, the result 
of the difference in SB profiles between Fl lOW and H interior 
to 2". 

4. DISCUSSION 

The new features we have observed in the disk at 1.1 /im 
add questions to the nature of the disk and the dust. The 
increase in PA towards the inner part of the disk could be 
caused by a genuine warp, or perhaps the presence of a sec- 
ond disk inclined at >12° from the main disk as is seen in 
/3-Pictoris. This is supported by the detection of an asymme- 
try in the northern part of the disk seen in the ACS and H data 
ofKFG07. 

KFG07 noted an overall blue V-H color for the Western 
lobe of the disk. The presence of a neutral scattering effi- 
ciency out to 1.1/tm as well as a red FllOW-H color inte- 
rior to 2" suggests that the spectrum of the HD 15115 disk is 
more complicated than the V and H data suggest. The abrupt 
change in scattering efficiency beyond 2" from FllOW to H 
may be due to the presence of strong absorption at ^1.65/im. 
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One candidate would be strong absorption due to water ice, 
something not yet detected in the scattered Hght of a disk. If 
that is the case, even deeper absorption should be present at 2 
and 3.5/im. Two possibilities exist to explain the red Fl lOW- 
H color of the inner disk-either there is evidence of strong 
absorption due to olivine (which has an absorption feature at 
1 /xm), or the presence of a very red materi al in the disk, possi- 
bly similar to what is seen in HP 100546 dArdila et al.ll2007h 
or HR 4796A (iDebes et al.ll2008l) . 

The age of HD 15115 is not well constrained, which 
would be helpful in understanding the origins of the debris 
disk. Kinematically, it is marginally consistent with mem- 
bership in the 12 M yr-old /3-Pictoris Moving Group (BPMG) 
dMoor et alj|2006l) . However, this is not borne out by back- 
tracking the position of HD 15115 based on its proper motion 
and radial velocity usi ng the methods of lOrtega et al.l (12002b 
and ISong et all (l2003jK_ nor is kinematics a reliable selector 
of a moving group Tsong etalJ[2002l). Most other indica- 
tors point to a > 100 Myr age. iNordstrom et al.l ( 120041) de- 
termined an age of 900 Myr, but this age is uncertain given 
their methodology-the lower limit for the age is 0, while the 
upper limit is 2.2 G yr. Ca II H and K line indicators point to 
an age of 500 Myr ( S il verston^ l2000l) . Although not an ac- 
curate age indicator for early F-type stars, the Lithium 6708A 
feature shows an equivalent width of 40 mA (Song 2008, in 
prep.), which is consistent with '^lOO Myr Pleiades stars. As 
a comparison, previously recognized 30 Myr old (Tucana- 
HorA) stars or younger early F-type stars all show equivalent 
widths for Li of ~ 100mA. The evidence to date points to an 
older age for HD 15115. 

HD 15115 is a prime example of the potential complexity 
of nearby circumstellar debris disks. It possesses a strange 
brightness asymmetry in its outer reaches that extends into 
about 2"(90 AU). This asymmetry, however, is a function of 
wavelength, and at least part of the cause may be composi- 
tional or grain population differences between the two sides 
of the disk. KFG07 investigated whether the close passage 
of a star might explain the strange morphology of the disk. 
They found a candidate object, HIP 12545, that might have 



passed very close to the debris disk of HD 15115, but while 
its motion is consistent with being a co-moving object, back- 
tracking of the the two star's positions back in time based 
on their Hipparcos/Tycho-2 astrometry and radial velocities 
shows that they are more widely separated in the past than 
they are now and could not have interacted. 

The fact that we measure a warp only in the western part of 
the disk around 1" suggests that HD 15115 has an even more 
complex structure. It would be useful to confirm the presence 
of this structure in the disk at other wavelengths. This is not 
easily done in the visible, unless STIS is refurbished in the 
next servicing mission-it is the only instrument able to image 
to within ^^1" in the visible. Follow-up H and K (or F160W or 
F205W images with HST) will be useful to look for the warp 
as well as constrain the ultimate nature of the composition of 
the grains in the disk. If water ice is indeed present, it will 
show strong absorption at around 2 and 3.5/im. Therefore, K 
and L band images of the scattered light disk are cruicial. 

What are the other possibilities to explain both the warp, 
the asymmetry, and the varying nature of the dust grains? One 
possibility is that the asymmetry is caused by a relatively re- 
cent collision that has injected many particles recently into the 
disk. Such an event is rare given the collisional times of the 
parent bodies and the timescale for such an event to be evenly 
spread in azimuth. Similar asymmetries have been seen in 
the mid-IR with /9-Pictoris ( Weinberger et al. 2003), and it is 
striking that HD 15115 shares many of traits with this well 
known disk. 
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